Coulomb breakup effects on the elastic cross section 
of ^He+^°^Bi scattering near Coulomb barrier energies 

T. Matsumoto,''QT. Egami,^ K. Ogata,^ Y. Iseri,^ M. Kamimura,^ and M. Yahiro^ 

'institute of Physical and Chemical Research (RIKEN), Hirosawa 2-1, Wako, Saitama 351-0198, Japan 
^Department of Physics, Kyushu University, Fukuoka 812-8581, Japan 
'^Department of Physics, Chiba-Keizai College, Inage, Chiba 263-0021, Japan 
(Dated: February 9, 2008) 

We accurately analyze the "He+^^^Bi scattering at 19 and 22.5 MeV near the Coulomb barrier energy, using 
the continuum-discretized coupled-channels method (CDCC) based on the n+n+'*He+^'^^Bi four-body model. 
The three-body breakup continuum of ''He is discretized by diagonalizing the internal Hamiltonian of ''He in 
a space spanned by the Gaussian basis functions. The calculated elastic and total reaction cross sections are in 
good agreement with the experimental data, while the CDCC calculation based on the di-neutron model of ''He, 
i.e., the n+ He+ Bi three-body model, does not reproduce the data. 
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In the recent measurements of ^HeH-^^^Bi scattering at 19 
and 22.5 MeV near the Coulomb barrier energy f?, '3L large 
enhancement of the a-emission cross section, which domi- 
nated the total reaction cross section, compared with that for 
the corresponding ^Li-induced reactions was reported. In or- 
der to clarify the nature of the enhancement, Keeley et al. 
analyzed the scattering by means of the continuum-discretized 
coupled-channels method (CDCC) that is a fully 

quantum-mechanical method for describing scattering of a 
three-body system. In the analysis, the ^HeH-^"^Bi system was 
assumed to be the ^nH-'*HeH-^°^Bi three-body system, that is, 
the neutron pair in ^He was treated as a single particle, di- 
neutron Cn). They found that the enhancement of the total 
reaction cross section of the ^HeH-^°^Bi scattering was due 
to the electric dipole {E\) excitation of ^He to its continuum 
states, i.e. Coulomb breakup processes of ^He, which was 
approximately absent in the ^LiH--^"^Bi scattering case. Their 
calculation, however, did not reproduce the angular distribu- 
tion of the measured elastic cross section and overestimated 
the measured total reaction cross section by a factor of three. 
Thus, reaction mechanisms of the ^HeH-^°^Bi scattering are 
not fully understood. 

In the very recent work fs*], it was reported that the elastic 
cross sections of the ^HeH--^"^Bi scattering calculated within 
the same framework as in Ref. |3] with the strength of 
the dipole coupling potentials multiplied by 0.5 reproduced 
the experimental data. This indicates that the E\ excita- 
tion strength of ^He cannot be accurately reproduced by 
the ^nH-'*He model. Since ^He is well known as a two- 
neutron halo nucleus, its structure should be described by 
the nH-n-H'^He three-body system rather than the ^nH-'*He two- 
body one. Thus, it is necessary to analyze the ^HeH-^^^Bi scat- 
tering by using the nH-nH-^HeH-^*'^Bi four-body model. Fur- 
thermore, a fully quantum-mechanical method such as CDCC 
is highly required to analyze the scattering near the Coulomb 
baiTier energies in which both nuclear and Coulomb breakup 
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processes are significant. 

In our previous work we proposed four-body CDCC 
that is the extension of CDCC and describes four-body 
breakup processes. In four-body CDCC, the three-body 
breakup continuum of the projectile is discretized by diago- 
nalizing the internal Hamiltonian in a space spanned by the 
Gaussian basis functions. So far the Gaussian basis func- 
tion was used with success for solving bound-state prob- 
lems of few-body systems. The approach is called the Gaus- 
sian expansion method (GEM) |10|. The application of 
the Gaussian basis function to the discretization of breakup 
continuum is a natural extension of GEM. In general, the 
method that describes the breakup continuum by a superpo- 
sition of i^-type basis functions is called the pseudostate 
method LI 1.1 . and other basis functions have been also pro- 
posed so far O d Q. Four-body CDCC was success- 
fully applied to ^He-n^^c scattering at 18 and 229.8 MeV in 
which only nuclear breakup was significant. The elastic and 
breakup cross sections calculated with four-body CDCC are 
found to converge as the number of Gaussian basis functions 
is increased 19]. This indicates that the set of discretized con- 
tinuum states obtained by the pseudostate method with GEM 
forms an complete set with good accuracy in a finite region 
of space that is important for the four-body reaction process 
concerned. Furthermore, in Ref. 1 15|, applicability of CDCC 
with the pseudostate method to Coulomb breakup processes, 
in which large modelspace of the projectile is required be- 
cause of the long-range property of Coulomb coupling po- 
tentials, was shown for the ^B-n'^^Ni scattering at 25.8 MeV. 
Thus, it is expected that four-body CDCC with GEM accu- 
rately describes the ^He-n^'^^Bi scattering near the Coulomb 
baiTier energy. 

In this Rapid Communication, we analyze ^He-n^'^^Bi scat- 
tering at 19 and 22.5 MeV by means of four-body CDCC. This 
is the first application of four-body CDCC to low-energy scat- 
tering in which both nuclear and Coulomb breakup processes 
are significant. We show that four-body CDCC reproduces 
the measured elastic and total reaction cross sections reason- 
ably well. Effects of four-body breakup processes on the elas- 
tic scattering are investigated through the dynamical polariza- 
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tion (DP) potential. We discuss the reason why the di-neutron 
model of ®He is insufficient to describe the ^He+^"^Bi scat- 
tering. 

We assume that the ^He+^*'^Bi scattering is described as the 
n+n+^He+^"^Bi four-body system. The model Hamiltonian 
of the system is defined by 

H = Kr + ?7„Bi(-Rrii) + C^nBil^na) + C^QBi(^a) + Hq,{\) 

where Rx (R) is the coordinate of particle X (the center- 
of-mass of ^He) relative to ^°^Bi, Kn is the kinetic energy 
associated with R, and Hq is the internal Hamiltonian of ^He. 
The potential C/„Bi (t^aBi) represents the interaction between 
n C'He) and •^"^Bi. It should be noted that UaBi contains a 
Coulomb part that causes Coulomb breakup processes in the 
^HeH-^°^Bi scattering. In four-body CDCC, the total wave 
function of the four-body system is expanded in terms of a 
finite number of the internal wave functions of the ^He pro- 
jectile. The internal wave functions including the bound and 
discretized-continuum states are generated with GEM as men- 
tioned above. In GEM, the n th eigenstate ^nim of ®He with 
the total spin / and its projection on the z-axis m is written as 
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where c denotes a set of Jacobi coordinates shown in Fig. [2 

(c) 

Each V'n/m expanded in terms of the Gaussian basis 

functions: 
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where A (£) is the angular momentum regarding the Jacobi 
coordinate yc (rc), and 77^/2 is the spin wave function of each 
valence neutron (n\ or 712). In actual calculation we truncate 
A and t at appropriate maximum values, A^ax and ^,„ax, re- 
spectively. The Gaussian range parameters are taken to lie in 
geometric progression: 



= fi(f,„ax/ri)(^-')/(-'''"^''-'^ 
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The eigenstate of ^He is antisymmetrized for the 

exchange between n\ and 712; we then have ^^^^"^5 — 

(-)^^aj"AS ™d (-)^^^ = 1 for c = 3. Meanwhile, the 
exchange between each valence neutron and each nucleon in 
''He is treated approximately by the orthogonality condition 
model 1 16]. The eigenenergies e„/ of ^He and the correspond- 
ing expansion-coefficients A|^j"as ^'^^ determined by diago- 
naHzing i/g 1 17ul8il. 

Using the internal states of ^He thus obtained, we expand 
the total wave function of the 7T,H-7H-'*HeH-^"^Bi four-body sys- 
tem with the total angular momentum J and its projection on 




FIG. 1 : Jacobi coordinates of the three rearrangement channels (c - 
1-3) adopted for the n+n+''He model of ®He. 
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jAf and L is the orbital an- 



gular momentum regarding R; below we denotes the channels 
{n,I,L} as 7. The expansion-coefficient represents the 
relative motion between the ^He projectile and the ^"^Bi tar- 
get and Pni is the corresponding relative momentum. Multi- 
plying the four-body Schrodinger equation {H — E)'^'^'^^ — 
by 3^*/*^ from the left and integrating over all variables ex- 
cept R, one obtains the set of coupled differential equations 
forx:! 
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where the coupling potential U.y'^{R) is defined by 

Uy^iR) = {y^,^'\UnBi{Rm) + UnBi{Rn,) + U^BiiRc.)\y^^') 

and /i is the reduced mass between ^He and ^"^Bi. We obtain 
the elastic and discrete breakup 5-matrix elements by solv- 
ing Eq. Q under an appropriate asymptotic boundary condi- 
tion |4, 19]. Details of the formalism of CDCC are shown in 
Ref. 14]. 



TABLE I: The maximum internal angular momenta and the Gaussian 
range parameters for each Jacobi coordinate. 
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3 0+ 
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0.1 


15.0 
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1,2 0+ 


1 
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15.0 
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15.0 


3 1- 
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0.1 


15.0 


0.5 


15.0 


1,2 1" 


1 
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0.5 


15.0 


0.5 


15.0 


3 2+ 


2 
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0.1 


15.0 


0.5 


15.0 


1,2 2+ 


1 


1 


0.5 


15.0 


0.5 


15.0 



In the present four-body CDCC calculation for ^HeH-^°^Bi 
scattering at 19 and 22.5 MeV, we take = 0+, 1", and 2+ 
states for ^He. Inclusion of 1^ state is essential to describe 
Coulomb breakup processes. As for the internal Hamiltonian 
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of ^He, we adopt the same Hamiltonian as used in Ref. j^l- 
We show in Table|I]the maximum values of the internal angu- 
lar momenta, Amax and ^max, and the Gaussian range param- 
eters, j/i, ijniax, ri, and fmax, used in the calculation of $„/m- 
For each set of {c,X,£,A,S}, we take imax = jms^x ~ 10. 
These values of the parameters are found to give good conver- 
gence of the calculated elastic and total reaction cross sec- 
tions. It should be noted that the maximum value of each 
Gaussian range parameter is 15 fm, which is quite larger than 
that used in the four-body CDCC analysis of ^He nuclear 
breakup i.e. 10 fm. Some parameters shown in TableU 
depend on I'^ and c, while in Eqs. Q-© the dependence has 
not been shown for simplicity. 

We select the $„/m with e„7 < 7 MeV among those ob- 
tained by diagonalizing and use them in actual CDCC cal- 
culation, since the ^nUn with e„/ > 7 MeV are found to have 
no effect on the calculated cross sections of the ^HeH-^°^Bi 
scattering shown below. The resulting number of the discrete 
states for the 0+, 1^, and 2+ states is 37 (including the ground 
state of ^He), 44, and 53, respectively. As for the nuclear parts 
of C/„Bi and UaBi, respectively, we take the optical potentials 
of Koning and Delaroche |20] and of Barnett and Lilley |2^. 
The maximum value of L is taken to be 200 and the scattering 
wave function is connected to its appropriate asymptotic 
form at i? = 200 fm. 

Below we show also results of CDCC calculation based 
on the ^n+'^He+^^^Bi three-body model to see three-body 
breakup effects on the elastic and total reaction cross sec- 
tions. In this model the di-neutron (^n) model of ^He is used 
and the intrinsic spin and the relative energy of '^n are as- 
sumed to be zero. We henceforth call CDCC based on the 
model above three-body CDCC. As for the interaction be- 
tween and ''He, we take the same parameter as used in 
Ref. Oil. The ^n-n^He continuum is discretized by the pseu- 
dostate method described in Ref. fllll and truncated at the rel- 
ative momentum k = 0.7 fm~^ that corresponds to about 7 
MeV of the excitation energy of ^He. The number of the 
discrete state is 9 for each of the 0"*", 1~, and 2+ states of 
^He. The range parameters of the Gaussian basis functions 
are (ai = 0.5 im,an== 20.0 fm,na = 20) with the same 
notation as in Ref. 111]. We adopt the optical potential of 
Barnett and Lilley |21] for the nuclear part of UaBi, as in the 
four-body CDCC calculation, and the d-'^^^Ph (type-a) optical 
potential at the deuteron incident energy of 15.0 MeV |22] for 
the interaction between and ^"^Bi. Other parameters of the 
modelspace are the same as in the four-body CDCC calcula- 
tion. 

Figure|2lshows the angular distribution of the elastic differ- 
ential cross section for the ^HeH-^*'^Bi scattering at 19 MeV. 
The solid line is the result of four-body CDCC and the dashed 
line is that of three-body CDCC. The results of four-body 
CDCC and three-body CDCC without breakup effects of ''He 
are shown by the dotted and dot-dashed lines, respectively. 
The difference between the solid and dotted (dashed and dot- 
dashed) lines shows effects of the four-body (three-body) 
breakup on the elastic cross section. One sees that four-body 
CDCC reproduces the experimental data well, while three- 
body CDCC underestimates the data at middle angles of 50°- 
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FIG. 2: Angular distribution of the elastic differential cross section as 
the ratio to the Rutherford cross section for the "He+^^^Bi scattering 
at 19 MeV. The solid (dashed) and dotted (dot-dashed) lines show the 
results of the four-body CDCC (three-body CDCC) calculation with 
and without breakup effects, respectively. The experimental data are 
taken from Ref. 1 1, 2]. The incident energy for the experimental data 
in the laboratory frame is shown to be 19 MeV and 18.4 MeV in the 
first ^lil and second \2J papers of Aguilera et al, respectively; in the 
present study we take 19 MeV. 
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FIG. 3: The same as in Fig.|2lbut for ''He+^°^Bi scattering at 22.5 
MeV. The experimental data are taken from Ref. fUjl^]. We take the 
incident energy of 22.5 MeV shown in the first paper of Aguilera et 
al. in. 

100°. The dashed line is consistent with the result of Ref. f3\; 
it should be noted that the real part of each coupling potential 
was multiplied by 0.8 in Ref. |3j], while in the present study 
such a renormalization factor is not included. Figure|3]shows 
the result at 22.5 MeV and features of the result are just the 
same as at 19 MeV. 

We show in Fig. |3 the calculated total reaction cross sec- 
tions; the solid (open) circles represent the results of four- 
body (three-body) CDCC. The open squares are the experi- 
mental values of the a-emission cross sections [IJ added by 
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FIG. 4: Total reaction cross sections of He+ Bi scattering as a 
function of the incident energy in the center of mass frame. The four 
values on the left (right) side correspond to the ®He incident energy 
of 19 MeV (22.5 MeV) in the laboratory frame. The solid and open 
circles represent the results of the four-body CDCC calculation and 
the three-body CDCC calculation, respectively. The experimental 
data shown by the open and solid squares are taken from Ref. 
see the text for details. 



the fusion cross sections |23|, while the soUd squares are the 
total reaction cross sections evaluated from an optical-model 
analysis of the measured elastic cross sections 1 1]. The figure 
shows that three-body CDCC overestimates the experimental 
data by a factor of about three. This overestimation is con- 
sistent with the fact that three-body CDCC underestimates the 
elastic cross section as shown in Figs. |2] and |3] On the con- 
trary, four-body CDCC reproduces the experimental data quite 
well, which shows the importance of the accurate descrip- 
tion of three-body breakup continuum of ^He. The remain- 
ing difference of about a few tens of % between the results 
of four-body CDCC and the data needs further investigation, 
including analysis of ^He-n^'^^Pb scattering in barrier-energy 
region |24]. 

In order to clarify the reason why the total reaction cross 
section calculated with three-body CDCC is much larger than 
that with four-body CDCC, we see first the strength of El 
transition based on the two models. The non-energy weighted 
El excitation strength B{E1) from the ground state $000 of 
^He to its excited states ^nim with / = 1 is given by 



B{El) 



^^\{^nlm\VsYlt,{yz)\^ 



000/ 



(8) 



The summation over n is taken up to a value that corresponds 
to the excitation energy of ^He of 7 MeV, i.e. the maximum 
energy of the modelspace of the present analysis. The re- 
sulting values of B{E\) based on the ^n-n^He model and the 
n+n+^Wt model are, respectively, 1.5 e^fm'^ and 0.9 e^ivc? . 
The latter agrees well with the experimental value reported by 
Aumann et al. |25], which is consistent with the conclusion 
in Ref. |26]. Thus, the El strength is overestimated in the 
di-neutron model. This indicates that the di-neutron model 



overshoots the breakup cross section of the ^He-n^^^Bi scat- 
tering, which is the main reason why three-body CDCC over- 
estimates the result of four-body CDCC, hence the measured 
total reaction cross section. This conclusion is qualitatively 
consistent with that drawn in Ref. |8], in which, as mentioned 
above, three-body CDCC with the strength of the dipole cou- 
pling potentials multiplied by 0.5 was shown to reproduce the 
elastic scattering data of the ^HeH-^°^Bi scattering. 

Next we discuss the difference of the optical poten- 
tials taken in the ^nH-^HeH-^"^Bi three-body model and the 
n+n+'^Wt+^^'^'Qi four-body model. As one sees in Figs.|2]and 
|3l the elastic cross section without breakup effects calculated 
with the three-body model underestimates that with the four- 
body model. This means that the diagonal component of the 
imaginary potential in the elastic channel based on the three- 
body model is deeper than that on the four-body model. Thus, 
the three-body model yields an absorption cross section larger 
than the four-body model does. This is also an important fac- 
tor for the enhancement of the total reaction cross section cal- 
culated with three-body CDCC. 

Finally, we see the difference between three-body CDCC 
and four-body CDCC in more detail by evaluating the dynam- 
ical polarization (DP) potential C/^p and the equivalent local 
potential U;!^{R). Explicit form of the two potentials is given 
by 



(9) 



and 



ui^{R)^uMR) + u-,,noiR), 



- ttJ 



where the subscript denotes the incident channel. 




FIG. 5: The imaginary parts of the dynamical polarization (DP) po- 
tential and the equivalent local potential for the ®He+^"^Bi scatter- 
ing at 22.5 MeV with J = 10. The solid and dotted lines, respec- 
tively, represent the results of the DP and equivalent local potential 
calculated with four-body CDCC. The dashed and dotted-dash lines 
correspond to those calculated with three-body CDCC based on the 
di-neutron model of "He. 



Figure |5] shows the imaginary parts and VK^p of 

C/qp calculated with four-body CDCC (solid line) and three- 
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body CDCC (dashed line), respectively, for the ^He+^osBi 
scattering at 22.5 MeV. Also shown are the results of the imag- 
inary parts of calculated with four-body and three-body 

CDCC, i.e. Weq'^^ and Weq^\ It should be noted that the fig- 
ure shows the potentials only in the peripheral region, i.e., 10 
fm < i? < 30 fm, where the scattering wave has nonnegligi- 
ble values. The results correspond to J = 10 at which the par- 
tial reaction cross section becomes maximum. It is confirmed 
that the DP potential hardly depends on J around J = 10 in 
the peripheral region shown in Fig.|5l It should be noted that 
C/jjp defined by Eq. (|9j, or U^^{B), is not a smooth function 
since it is divided by the oscillating function x^g ■ Each Une in 
Fig.|5lis obtained by interpolation of [/j3p(i?) (or Ugq{R)) at 
R where \x:^„ \ is more than 90% of its maximum value in the 
asymptotic region; this simple way of evaluating t/^p (i?) and 
Ucq{R) is appropriate for the present purpose. The 5-matrix 
element calculated with the interpolated DP potential repro- 
duces the one obtained by CDCC within the relative error of 
3%. 

One sees from Fig. |5] that both W^Qp'^^ and W^^p'''' have a 
long ranged tail, which is induced by Coulomb breakup pro- 
cesses. In the tail region of i? > 15 fm, where the imagi- 
nai-y part of U^^jg is negUgible and Wfjp agrees with W^'^, 

W^p^ is deeper than W^p^'' . This is consistent with the fact 
mentioned above that the strength of El transition is overes- 
timated by the di-neutron model of ^He, which is used in the 

J(3) 

three-body CDCC calculation. On the other hand, W^p is 

shallower than Wj^p^^ in the region of i? < 15 fm. How- 
ever, the imaginary part of U^g^g calculated with three-body 
CDCC is, as mentioned above, much deeper than that with 

four-body CDCC, which makes Wcq^"^ deeper than Wcq'^^ in 
this region. Thus, the imaginary part of the equivalent lo- 
cal potential, which dictates the total reaction cross section, 
calculated with three-body CDCC is deeper than that with 
four-body CDCC in the entire region that is important for the 
^HeH-^°^Bi scattering concerned. This conclusion is consis- 
tent with the discussion above and the result shown in Fig.0] 
The real part V^p of the DP potential also plays impor- 
tant roles in the ^HeH-^"^Bi scattering. It is found that Vp(p 



is repulsive for i? < 15 fm, where nuclear breakup processes 
are significant, and attractive for i? > 15 fm, where Coulomb 
breakup processes are dominant. The difference between V^p 
calculated with three-body CDCC and four-body CDCC in the 
outer region is just the same as for Wp,p shown in Fig.|5l 

In summary, the ^He-H^^^Bi scattering at 19 MeV and 
22.5 MeV, near the Coulomb barrier energy, is analyzed with 
the continuum-discretized coupled-channels method (CDCC) 
based on the nH-rt-n "'Hen- ■^"^Bi four-body model, four-body 
CDCC, that treats both nuclear and Coulomb breakup pro- 
cesses simultaneously. It is found that the angular distribu- 
tion of the elastic cross section and the total reaction cross 
section calculated with four-body CDCC reasonably repro- 
duce the experimental data with no free adjustable parame- 
ter Three-body CDCC based on the '^ii+'^He+^^^Bi three- 
body model, i.e. di-neutron model of ^He, turns out to over- 
estimate the total reaction cross section by a factor of about 
three, which is consistent with the conclusion of the previ- 
ous work of Keeley et al. |,3]. The value of B{E1) calcu- 
lated by the di-neutron model overestimates the experimen- 
tal one, while the n+n+'^He three-body model reproduces the 
data well, as pointed out in Ref. 1 8]. This makes the imaginary 
part of the dynamical polarization potential calculated with 
three-body CDCC deeper than that with four-body CDCC in 
the tail region. Also important fact is that three-body CDCC 
contains a deep imaginary part of the diagonal potential for 
the elastic channel compared with four-body CDCC does. We 
thus conclude that the di-neutron model is inadequate to treat 
breakup continuum of ^He precisely and a n+n+'^He+^'^^Bi 
four-body reaction model is necessary to accurately describe 
the ^HeH-^*'^Bi scattering. Four-body CDCC is indispensable 
to analyze low energy ^He scattering in which both nuclear 
and Coulomb breakup processes are significant. 
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